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Abstract 

We consider the network model of the integer quantum Hall effect transi- 
tion. By generalizing the real-space renormalization group procedure for the 
classical percolation to the case of quantum percolation, we derive a closed 
renormalization group (RG) equation for the universal distribution of conduc- 
tance of the quantum Hall sample at the transition. We find an approximate 
solution of the RG equation and use it to calculate the critical exponent of 
the localization length and the central moments of the conductance distribu- 
tion. The results obtained are compared with the results of recent numerical 
simulations. 
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I. INTRODUCTION 



The fact that at zero temperature the quantum Hall transitions are infinitely narrow com- 
plicates strongly their theoretical description even without electron-electron interactions. By 
now, in theory, the quantitaive information about the localization of two-dimensional elec- 
tronic states in a magnetic field was obtained exclusively from numerical simulations. These 
simulations are carried out in one of two limits: 

i) Short-range disorder 

The results obtained before 1995 are reviewed in Refs. ^ and ^. Different approaches 
employed: the recursive Green's function method!, evolution of the number of states with 
non-zero Chern numbers with the sample size0, and finite-size scaling analysis of the Thou- 
less numberili give consistent results for the critical exponent, u, of the localization length. 
Direct calculations of the conductance at the transitionii show that it is close to 0.5e^/h. Re- 
cent developments include the study of the role of spin-orbit scattering for spin-degenerate 
electrons, localization in doublell- and multilayer systems!, statistics of energy levels at the 
transition^ and fractal properties of eigenstates at higher Landau levelJiil. 

a) Smooth disorder 

All the simulations^^ were performed in the frame of the network model introduced 
by Chalker and CoddingtonO. The questions addressed were essentially the same as in the 
case of the short-range disorder. In addition calculating the network model was used 

to investigate the spin-degenerate case0'0, and the effect of Landau level mixing on the 
positions of delocalized stateJlllll. It was also used to demonstrate that this mixing does 
not change the universality class of the transitionEl, to relate the longitudinal and transverse 
conductivities in the transition regionEl and to trace the smearing of a single delocalized state 
into a metallic band in a superlatticeS. 

Recent worksil^ explored the correspondence between the network model and certain 
limits of the models that were already studied (spin chain&i, Hubbard chainS and Dirac 
fermions^IEa) . 
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In the recent experimental paper by Cobden and Koga the observation of strong and 
reproducible conductance fluctuations of a mesoscopic sample right at the transition was 
reported. This observation was also accounted for in the frame of the network model in Refs. 
pO| and The result of simulations0 closely resemble the histogram of the conductance 
fluctuations inB In Ref. the calculated central moments of the conductance distribution 
were extrapolated to large sample sizes, where it becomes universal. 

The network model0 is illustrated in Fig. 1. If a disorder is smooth, classical electron 
drifts along equipotentials, which are closed and disconnected except for the percolation 
threshold. For simplicity it was assumed that at any energy equipotentials are arranged 
periodically. Quantum mechanics is introduced by allowing a strong tunnel coupling between 
the orbits when they come close enough (~ or closer than the magnetic length). This destroys 
the classical localization and makes important the interference of different directed paths. 
Randomness is introduced by assuming the magnetic phases, acquired in course of drift 
along equipotentials, to be random. 

If a smooth potential is symmetric with respect to zero energy, then the classical perco- 
lation threshold is also at e = 0. Although this is obvious, we would like to derive the same 
conclusion following a different line of arguments. 

The centers of closed drift trajectories in Fig. 1 form a square lattice. Classical electron 
switches to another trajectory, centered at a neighboring lattice site, if its energy, e, exceeds 
the height, Vq, of a saddle-point separating these sites. One can say that in this case a bond 
between two sites connects them. In the opposite case, e < Vq, the electron retains its orbit 
after passing the saddle-point. Then one can say that the corresponding bond is removed. 
Thus we arrive at the bond percolation problem on the square lattice. 

Suppose now, that the electron energy, which was initially at the lowest minimum, grad- 
ually increases. As it happens, more and more bonds switch on. Percolation occurs when 
some critical portion of bonds connect. It is an exact result of the percolation theoryil that 
for the square lattice this portion is 1/2, which corresponds to zero energy. 

The fact that the classical picture can be reformulated as a percolation problem on a 
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lattice suggests the following possibilty. It is known that the critical exponent for classical 
percolation can be found with high accuracy from a simple renormalization group (RG) 
calculation. In course of this calculation only a small cluster is considered. If a similar 
approach could be devised for the quantum percolation, it would provide an approximate 
description of the quantum Hall transition in the closed form. In contrast to the classical 
renormalization scheme, such an approach should deal not with "connecting" or "removed" 
bonds, but characterize each bond by some scattering matrix. It should also allow for 
interference of different directed paths within a cluster. This program is carried out in 
the present paper. We utilize a real-space RG scheme for classical percolation, proposed 
by Reynolds, Stanley, and Klein§, and generalize it to the quantum percolation. This 
generalization is described in Sec. II. We derive the RG transformation, which governs 
the evolution of the distribution of conductance right at the transition, with increasing the 
sample size. In Sec. Ill the approximate numerical solution for the fixed point of the RG 
transformation is obtained. This solution determines the universal conductance distribution 
at the transition. It is used for calculation of the exponent v in Sec. IV. In this section 
we also discuss the finite size corrections. In Sec. V our results for the distribution of 
conductance of a large sample at the transition and in its vicinity are compared with recent 
numerical simulations. Concluding remarks are made in Sec. VI. 



II. RENORMALIZATION PROCEDURE 



Let us briefly remind the renormalization procedure for the classical bond percolation 
problem on the square lattice!!^. In this procedure a fragment of original lattice containing 
8 bonds (Fig. 2a) is replaced by a fragment consisting of two superbonds. The probability, 
p, that a superbond does connect two supersites is expressed through the corresponding 
probability, p, for original lattice in a following way. Firstly, the supersites are connected 
if all five bonds, which constitute the superbond (see Fig. 2a), connect. The probability 
of this is p^. Secondly, removing one of five bonds, does not affect the connectivity of 
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the superbond, the total probability of such configurations being 5p^{l — p). Adding the 
probabilities for the superbond to connect when two and three original bonds are removed, 
one gets the following expression for renormalized probability. 

p=p^ + 5p^(l -p) + 8p^(l - pY + 2/(1 - pf. (1) 

The fixed point of the transformation (|1|), which determines the percolation threshold, is 
p = p = 1/2 and it coincides with the exact result. The equation for the critical exponent, 
i>c, of the correlation length, C,c, emerges as a condition that C,c is the same for original 
and for renormalized lattices. If the lattice constant of the original lattice is a, then = 
a(| —p)"'^"- On the other hand, the lattice constant of the renormalized lattice is 2a, so 
that = 2a(^ — p)^'^"- Assuming that — p) and — p) are small, the two relations 



yield u 



In 2/ \n{dp/dp) 



. Using Eq. (|T]), one gets Uc ~ 1.428, which is only 8 percent 

p=l/2 

bigger than the value 4/3, which is presumed to be exact. 

As it was mentioned in Introduction, we are interested in the situation when each bond 
represents a saddle-point of the random potential. Then we come to Fig. 2b, which can 
be viewed as a transformation of five original saddle-points into a super-saddle-point. In 
classical mechanics, p is the probability that the energy of an electron exceeds the height 
of the saddle-point. In quantum mechanics, a saddle-point is characterized by reflection 
coefficient, r, which is the amplitude for an incoming wave to retain the equipotential, or by 
transmission coefficient t-the amplitude to switch the equipotential. Obviously, + = 1. 
In the network, there are two incoming and two outgoing waves at each saddle-point i. The 
S-matrix, relating the amplitudes of outgoing waves to the amplitudes of incoming waves, 
has the form 



(2) 



With five saddle-points we have five matrix equations. Solving the system, we get the 
following expression for the effective transmission coefficient of the super-saddle-point 
tit5(r2r3r4e^^^ - 1) + hUe'^^^+'^^Knr^r^e-'^- - 1) + hjhhe'^^ + ^1^46^^^) 
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where $j are the random phases acquired by an electron after traversing one of four loops 
shown in Fig. 2b. This expression can be viewed as a generalization of Eq. ([1|) to the case 
of quantum percolation. It appears very convenient to parametrize ti, as follows 

- 1 _ 1 

'^^ - (e-. + 1)1/2' ''i - (e-.. + 1)1/2- (4) 

Then the parameter Zi can be related to the height of the saddle-point, Vi, asS 

F -V- 1^/2/2 
Zi = — — , r = — . (5) 

Here e is the energy of the electron, V^^, Vyy are the second derivatives of the saddle-point 
potential, and / is the magnetic length. Thus, for e = 0, Zi is the dimensionless height of 
the saddle-point and Eq. @ determines the dimensionless height, Z{zi,^j}, of the super- 
saddle-point for a given set of Zi and ^j. We now introduce the distribution function 

Kiz,{z,}) = {5iz-Z{z,,^,}))^^^y (6) 

of the heights of super-saddle-points at given heights of original saddle-points. The averag- 
ing in (|^) is carried out over the random phases The crucial test for applicability of the 
renormalization procedure to the quantum percolation is that the position of the delocalized 
state should remain unchanged after renormalization. Suppose that in original lattice all Zi 
were zero. Then the delocalized state is at e = 0. After the first step of renormalization the 
distribution of the heights of super-saddle-points is given by K{z, {zi = 0}). This function 
is plotted in Fig. 3a. We see that it is symmetric with respect to 2; = 0, which justifies the 
further renormalization. Note that Zi = at the first step, implies that the power transmis- 
sion coefficient, t^, is equal to 1/2 for all saddle-points. As it follows from Fig. 3a, after 
the first step we already get a substantial spread in (the characteristic range is 0.2+0.8). 
This is the result of interference of different paths, which is described by Eq. (^. 

Now we can write down the renormalization group equation. Denote with Qn{z) the 
normalized distribution function of the heights of super-saddle-points after the n-th step of 
renormalization. Then the coresponding distribution function after the step n + 1 is given 
by 



Qn+l{z)=f{Qr,}= /f n dzAK{z,{z,}) H Qni^i)- (7) 

i=l,..,5 i=l,..,5 

The fixed point, Q{z), of tlie transformation determines the distribution of conductance 
of a large sample at the plateau transition. The equation for the fixed point, Q{z) = T{Q}, 
is a quite complicated integral equation. In the next section we obtain its approximate 
solution. 



III. APPROXIMATE SOLUTION OF THE RG EQUATION 

The crucial observation which allows to obtain an approximate solution for the fixed- 
point distribution Q{z) is that within the relevant region of change of the variables of 
integration Zi, which appears to be J2 zf < 10, the shape of the function K{z^ {zi}) changes 
rather weakly. At the same time the position of the symmetry axis of K{z, {zi}) depends 
strongly on the values of Zi. In other words, the kernel of the transformation can be 
presented as K{z — f{zi\), where the function K{z) is calculated for all Zi = and shown 
in Fig. 3. By plotting K{z, {zi}) for many sets {zi} (some examples are shown in Fig. 3b) 
we had established an approximate form of the function /: 

f{Zi} = Ci{Zi + Z2 + Z4 + Zr,) + C3Z3 - d{zi - Za){z2 - Z^) - 

-\i{zl + zl + zl + zl) - X^zl + \2Z:i{zl + zl + zl + zl), (8) 

with coefficients ci ^ 0.416, C3 ^ 0.164, d ^ 0.067, Ai ^ 0.009, A2 ^ 0.009 and A3 ^ 0.003. 
The quadratic term in the expansion is non-physical. It results from the deficiency 
of the renormalization scheme, which is the absence of complete tx/2 rotational symmetry 
(see Fig. 2b). Fortunately, the corresponding coefficient, is small enough, so that this 
term has a small effect on Q{z). We have verifyed it by taking this term into account as a 
perturbation. This term also does not affect the value of the critical exponent, since it 
does not change upon a shift of all Zi by the same value (see below). We have also established 
(by first neglecting and then taking into account perturbatively) that the cubic term with 
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coefficient A2 leads to a relatively small correction to Q{z). On the contrary, the cubic term 
with coefficient Ai appears to be very important. 

After neglecting the quadratic and the last cubic terms in the function f{zi} decom- 
poses into the sum of functions of different arguments. Owing to this, a great simplification 
can be achieved by taking the Fourier transform of both parts in Eq. (|^. One gets 



where Qni^^) is the Fourier transform of Qn{z) and K{uj) is the Fourier transform of K{z). 
Now the solution for the fixed point Qn{z) = Qn+i{z) = Q{z) can be easily obtained 
numerically. The result is shown in Fig. 4. Except for the tails at |z| > 5 it can be well 
approximated by a gaussian ^J^exp{—az'^) with a ~ 0.1. In the right-hand side of Eq. (|^ 
there is a product of five functions Q{zi). Having a ~ 0.1, the 5-dimensional volume, which 
provides the major contribution to the integral in (^, can be estimated as J2 zf < 10. We see 
that the fixed-point distribution is broad. The half-width is z ~ ±2.5, which corresponds 
to the range 0.08 4- 0.92 for the power transmission coefficient. In the next section we use 
the solution obtained for calculation of the critical exponent u. 



IV. CRITICAL EXPONENT 

We first address the critical exponent of the localization length. Suppose that before the 
renormalization the distribution of z is already given by Q{z). Then for electron with a small 
but non-zero energy, e, this distribution is given by Q{z — zq), where zq = e/T <^ 1 (see Eq. 
(D). Substituting Q{z — zq) into the right-hand side of the transformation one gets a 
renormalized distribution, Q{z — tzq), where r is some number, independent of zq. Upon 
repeating the procedure, the shift will grow as zqt^ and, eventually, become ~ 1. At this 
step an electron will be strongly localized in a renormalized lattice of super-saddle-points 
(within ~ 1 unit cell). Thus, one should identify the localization length, ^, with a lattice 
constant after these n steps, which is equal to 2"a. Then the condition -Zot" ~ 1 can be 
rewritten as ^(i)^''^'^/^"^) 1^ and one has ^ ~ a(^)~'^ with 
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In2 , , 

The remaining task is to calculate parameter r using Eq. (^. For small Zq the result of 
substituting Q{z — zq) into the right-hand side can be presented as 

Qi^) + ^0 /( n d^^) (e ^) n Qi^^) (11) 

\=1..5 ^ ^ j ^ i=1..5 

At small z the first term behaves as Q{0) + yQzz{0), while the second term is proportional 
to z (due to symmetry). As a result, the maximum of the sum ( pT]) will get shifted from 
z = to z = ZqT, where 

1 /■/ „ . \ d^K 

T = 

To proceed further we adopt the same simplification, K{z,{zi}) = K{z — f{zi}), as was 
adopted to find the fixed point Q{z). Then one has 



J2 TTir- 1^=0= irri-fi^i}) 4ci + cg - 3Ai(zi2 + + + z^) - SX^z^ 



(13) 



Note that non-physical quadratic term in (§) did not contribute to (p^ . After substituting 



(|13|) into ( [T2I) it is again convenient perform the Fourier transformation and then to make 
use of Eq. (^ . Finally we get 

^ = ^'^ + '^"^ Jdu;u^Q{u) ' 



where 



Evaluating numerically the integrals in ( plSj) and then in (p!^ with (^(-z) found in the previous 
section, we got for r the value r = 1.336, which leads to the exponent z/ = 2.39. Since Ai, A3 
were determined with certain error bars, we repeated the calculations and for Ai = 0.01 (the 
effect of A3 is small) and got r = 1.306, u = 2.59. The results of most numerical simulations 
are grouped around z/ = 2.4. 



Despite the agreement obtained, we cannot estimate the accuracy of the result, since the 
calculation was based on the assumption that the shape of the kernel of the transformation 
(0) is independent of Zi. In fact the width of the distribution K{z, {zi}) decreases slowly 
with increasing \ Zi\, and, correspondingly, Kiuj) broadens. We found that at the boundary 
of the relevant region of the change of Zi, ^zf < 10, the width of K{uj) increases by ~ 30 
percent. We could not include this effect in our calculation. To obtain the limiting estimate 
we repeated the calculation assuming that the width of K{uj) is 30 percent bigger for all Zi. 
This caused the reduction of z/ to 1.95. 

In conclusion of the section let us discuss the question how the distribution Qn{z) ap- 
proaches the fixed-point distribution Q{z) with increasing the sample size. To answer this 
question let us introduce the deviation SQn{z) = Qn{z) — Q{z) and linearize the transfor- 
mation (0). After using the above simplification for K{z, {zi}) and performing the Fourier 
transform, we get 



5Qn+i{co) = F{5Qn} = K{uj)Jl{uj) J dz 



(16) 

where 

j,_3(o;) = / dzQ{z)e'''^''^''-^'''''\ (17) 



If initial perturbation was 6Qq{z), then after n steps it will evolve to F"{6Qq{z)}. The 
relevant perturbations are the even functions of z, satisfying the condition / dzQQ{z) = 
(i. e. Qo{uj ^0) = 0), to keep Qn{z) normalized. We were able to make the complete 
analysis of the transformation (0) only in the case Ai = A3 = 0. The answer is as follows. 
The fluctuations for which, at — > 0, 6Qo{uj) behaves as dj^™, m being integer, decay with 
increasing n as {Acf^ + cf^)"^ oc L"'^™, where 7^ = In^Acf^ + cf^)'^ / \n2. The slowest decay 
is characterized by the exponent 71 ~ 0.5. 

With finite Ai,A3 there is no such a simple classification. We could only see that the 
exponent increases significantly up to the value ~ 2. At large L this exponent would 
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determine the finite-size correction to the conductance distribution and, thus, to the average 
conductance, variance and all the moments. 



V. CONDUCTANCE FLUCTUATIONS 

The fixed-point distribution Q{z) is directly related to the distribution function of the 
conductance at the plateau transition. After an appropriate number of the renormalization 
steps the entire sample reduces to a single super-saddle-point; the power transmission co- 

2 

efficient of this super-saddle-point would determine the conductance of the sample, ^G, 
with G = (e^ + 1)^^. Then the distribution function of G can be found from the relation 
P{G)dG = Q{z)dz yielding 

The distribution function P{G) is shown in Fig. 5. It is very broad. This fact was first 
pointed out in the experimental paper^l and then reproduced in the numerical simulations 
of Cho and Fisher§. As a distinctive feature, the distribution function (^) has a shallow 



minimum at G = 1/2. Using this function we calculated the variance 6G = y {G"^) — | and 
obtained the value 6G = 0.33, which agrees with the result of Ref. 0. To compare our 
results with the recent simulations by Wang, Jovanovic and D.-H. Led! we have calculated 
the moments of the distribution P{G), defined as An = {{G — |)"), where n is an even 
integer which we have changed from n = 2 up to = 20. The results are shown in Fig. 



6 together with the values A2,A4,Aq and As, obtained in the simulations of Ref. |3y. The 
agreement is quite good. 

Finally we have studied how the variance, 6G{e), falls off as the energy, e, deviates from 
the position of the delocalized state e = 0. As it was mentioned above, for non-zero energy 
the distribution of z is given by Q{z — -^), where Eq = ^ , L being the sample size. 
Then one has 
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The dependence SG{e) is shown in Fig. 7. To compare it with the simulations of Ref. 



31] we had to estabhsh the correspondence between our dimensionless energy e/Eq and the 



one adopted in Ref. 0. We did it by calculating the average {G{e)) and fitting it to the 
curve in0. The agreement is again quite reasonable, especially taking into account some 
asymmetry of the numerical results. However, strictly speaking, the distribution of z at 
finite e has the form Q{z — -^) only if e/Eq is smaller than 1. 

Note in conclusion of the section, that identifying the conductance G with the transmis- 
sion through the last super-saddle-point implies that sweeping the electron energy through 
e = causes an increase of G from to 1. This was the case for the experimental geometry 
of Ref. m. 



VI. CONCLUSION 

Certainly, the renormalization scheme for quantum percolation, developed in the present 
paper, provides only an approximate description of the plateau transition. As in classical 
percolation, when, in the course of renormalization, one ascribes a certain probability to a 
superbond with no care about the condition of surrounding bonds, our quantum generaliza- 
tion leaves out many interference processes. This happens at the stage when we average over 
the phases, in Eq. (P). By doing so we assume the phases of transmission coefficients 
of super-saddle-points to be uncorrelated, which in reality is not the case. However, to the 
best of our knowledge, the above approach is the first one in which the characteristics of 
the quantum Hall transition are obtained from the solution of a certain closed equation. 
By taking into account only basic interference processes at each scale, we were able to re- 
produce the results of numerical simulations in which all interference processes were taken 
into account. Our approach does not require extrapolation to infinite sizes. This might be 
important in the cases when the result of such an extrapolation is rather ambigous. For 
example, when there are two (several) delocalized states very close in energy (Zeeman split 
levels, higher Landau levels). In this case the network should carry more than one chan- 
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nel and allow for their mixingEllllllilll0. Formal generalization of the RG equation to 
this case is straightforward, however its approximate solution (as in Sec. Ill) might be not 
sufficient. Note finally, that the RG results for the classical bond percolation problem can 
be significantly improved by considering a bigger clusterS'i. Using the rules, formulated in 
Sec. II, this advanced scheme can be also generalized to quantum percolation. 
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FIGURES 

FIG. 1. Schematic representation of the network model. Circles show equipotentials. Dashed 

lines show the saddle-points. 

FIG. 2. Schematic illustration of the renormalization procedure for the classical (a) and quan- 
tum (b) percolation. 

FIG. 3. Distribution of heights of the super-saddle-points K{z, {zi}) is plotted for dif- 
ferent sets {zi,Z2,Z3,Z4,Z5}; (a){0,0,0,0,0}; (b){-2,2,0,-l,l} (solid curve), {-1.5,-1.5,-2,-1.5,-1.5} 
(dotted curve), {-1.5,1.5,0,-1.5,-1.5} (dashed-dotted curve), {1,-1,2,1,5,0.5} (short-dashed curve) 
{1,2,2,1,2} (long-dashed curve). 

FIG. 4. Fixed-point distribution of heights of the super-saddle-points. 

FIG. 5. Distribution function of conductance at the quantum Hall transition. 

FIG. 6. Moments An = {{G — ^)") of the conductance distribution at the transition are shown 
for different n. Diamonds show the numerical results of Ref. 30. 

FIG. 7. Variance of the conductance fluctuations is plotted v.s. dimensionless energy. 
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